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Gravitational waves in cylindrically symmetric Einstein 
gravity are described by an effective energy tensor with the 
same form as that of a massless Klein-Gordon field, in terms 
of a gravitational potential generalizing the Newtonian poten- 
tial. Energy-momentum vectors for the gravitational waves 
and matter are defined with respect to a canonical flow of 
time. The combined energy-momentum is covariantly con- 
served, the corresponding charge being the modified Thorne 
energy. Energy conservation is formulated as the first law ex- 
pressing the gradient of the energy as work and energy-supply 
terms, including the energy flux of the gravitational waves. 
Projecting this equation along a trapping horizon yields a 
first law of black-hole dynamics containing the expected term 
involving area and surface gravity, where the dynamic surface 
gravity is defined with respect to the canonical flow of time. A 
first law for dynamic cosmic strings also follows. The Einstein 
equation is written as three wave equations plus the first law, 
each with sources determined by the combined energy tensor 
of the matter and gravitational waves. 
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I. INTRODUCTION 

The emission of gravitational waves by black holes is 
under intense investigation, in preparation for the ex- 
pected observational study of black holes and other as- 
trophysical objects by upcoming gravitational-wave de- 
tectors. Most of this work is either numerical or in an ap- 
proximation, since we lack a full physical understanding 
of the dynamics of black holes and gravitational waves. 
The accepted theory of black holes mainly concerns stat- 
ics and asymptotics; for instance, there are standard 
asymptotic definitions of mass, angular momentum and 
the energy flux of gravitational waves, with surface grav- 
ity defined for stationary black holes. One might expect 
such physical quantities to play a more local role in dy- 
namical processes. 

A framework for addressing black-hole dynamics was 
introduced a few years ago Q. Black holes were given 
a local, dynamical definition in terms of trapping hori- 
zons, which are essentially the locations where light waves 
are marginally trapped by the gravitational field. Some 
basic properties of black holes were established, such as 
that the horizon is achronal and therefore locally one-way 
traversible, and a second law expressing the increase of 
area of the trapping horizon. A comprehensive picture of 



the spherically symmetric case has since been developed, 
involving local, dynamical definitions of physical quanti- 
ties such as energy Q and surface gravity || , related by 
local, dynamical equations including a first law. 

This article introduces the corresponding physical, ge- 
ometrical quantities and equations in cylindrical symme- 
try, which has the additional complexity of gravitational 
waves. In vacuo, these are known as Einstein-Rosen grav- 
itational waves [Q. Of course, cylindrical black holes are 
rather unphysical, but their study is a precursor to that 
of axially symmetric black holes, in particular their in- 
teraction with gravitational waves. The main conceptual 
issue is whether the energy flux of gravitational waves 
admits a local, dynamical definition. This turns out to 
be so, indeed quite natural, in that there is a covariant 
conservation law for the combined energy-momentum of 
the gravitational waves and matter. 

Due to some peculiar properties of black holes in cylin- 
drical symmetry, naked singularities may also be a fea- 
ture of gravitational collapse. Consequently the frame- 
work is extended to include axial singularities, particu- 
larly dynamic cosmic strings. 

The article is organised as follows. Section II reviews 
cylindrical symmetry in the sense of Melvin |g and de- 
fines regular axes and axial singularities. Section III de- 
fines black holes and related ideas of gravitational trap- 
ping. Section IV reviews the modified Thorne energy 
[Q, which turns out to be the appropriate definition of 
gravitational energy. Section V introduces a canonical 
flow of time and the corresponding energy-momentum 
vector of the matter. Section VI introduces the energy- 
momentum vector of the gravitational waves, states the 
conservation law, and introduces the gravitational po- 
tential and effective energy tensor of the gravitational 
waves. Section VII states the first law in terms of the 
energy flux of the gravitational waves, plus energy-flux 
and work terms for the matter. Section VIII defines the 
surface gravity. Section IX collects basic laws of black- 
hole dynamics, including the projection of the first law 
along a trapping horizon, which involves area and surface 
gravity as expected. Section X discusses some static ex- 
amples including black holes and cosmic strings. Section 
XI defines dynamic cosmic strings and discusses issues 
of cosmic censorship and predictability. Section XII con- 
cludes. An Appendix derives the Einstein equation and 
lists various expressions in standard coordinate systems. 
Einstein gravity is assumed, though the application of 
the Einstein equation is stated as appropriate. 
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II. REGULAR AND SINGULAR AXES 

There is sometimes discussion about how to distinguish 
cylindrical, planar and toroidal symmetry. Each may be 
defined locally by the existence of two commuting, spa- 
tial, Killing vectors, such that the orthogonal space is 
integrable. Then there exist coordinates (ip, z) such that 
the Killing vectors are (£ v ,£z) = (d/dip,d/dz). Here it 
will also be demanded that each of the Killing vectors 
be hypersurface orthogonal, which in planar symmetry 
is often called the polarized case. Then central cylin- 
drical symmetry may be defined by the existence of an 
axis, defined shortly, with respect to exactly one of the 
Killing coordinates, the azimuthal angle ip, which will be 
identified at and 2-7T. Central toroidal symmetry may 
be defined by the existence of axes with respect to both 
coordinates, and planar symmetry by no such axes. For 
local cylindrical, toroidal and planar symmetry, one re- 
quires only the identifications, without prejudice as to 
the existence of axes. This may be useful since, for in- 
stance, there may be black holes which look cylindrically 
symmetric from the outside, but have no axis inside. 

The norms of the Killing vectors are geometrical in- 
variants, the circumferential radius 

P = yV 5 (!) 
and the specific length 

t = (2) 

where the sign convention is that spatial metrics are pos- 
itive definite, the dot denotes contraction and the flat b 
denotes the covariant dual with respect to the space-time 
metric, i.e. index lowering. Similarly, a sharp jj will de- 
note the contravariant dual. Here and throughout the 
article, specific refers to the need to scale quantities by 
some specific length of the infinite cylinder, which is arbi- 
trary given the freedom to rescale £ z by a constant factor. 
Likewise, a cylinder will always mean a cylinder of sym- 
metry. Coordinates for the 2-dimensional quotient space, 
orthogonal to the cylinders, will not be taken in the main 
text, in order to stress the geometrical invariance of the 
quantities and equations. See the Appendix for a list of 
the required quantities expressed in a standard coordi- 
nate system. 

An axis will be defined as a boundary of the quotient 
space-time coinciding with p = 0, such that £r l = 0(1). 
Then £ is either finite and non-zero at the axis, or be- 
comes infinite. An axis will be said to be regular if £ is 
finite and 

Vp- V*p= 1 + 0(P 2 ) (3) 
Vp ■ = 0{p) (4) 
W-V<^ = 0(1) (5) 

and singular otherwise, where V is the covariant deriva- 
tive of the space-time metric. A singular axis will also 



be called an axial singularity. The first of the above con- 
ditions is standard, e.g. ||, but the others appear to be 
new. Some such conditions are necessary, since the first 
condition alone allows curvature singularities at the axis, 
e.g. for £ = 1 + p + 0(p 2 ). 

This article will adopt the above definition of central 
cylindrical symmetry for definiteness, though local cylin- 
drical symmetry suffices for some of the results. This is 
more rigid than a recent suggestion of Carot et al. || 
which does not require the orthogonal space to be inte- 
grable. One might call this twisted cylindrical symme- 
try. This is analogous to the relaxation from static to 
stationary space-times and would be expected to allow 
angular momentum. The present definition is essentially 
the whole-cylinder symmetry of Melvin |5|||, though it 
should be stressed that the axis has not been required to 
be regular. In particular, cosmic strings are allowed. 

III. BLACK HOLES 

The function 

r = p£ (6) 

plays a similar role to the areal radius of spherical sym- 
metry, with £ scaling various physical quantities so that 
they have the physically correct units. In particular, 

A = 2nr (7) 

is the specific area of the cylinders and 

V = vrr 2 (8) 

is the corresponding volume. The key ideas of gravita- 
tional trapping jl| may then be defined directly in terms 
of r or A, as in spherical symmetry ||, as follows. Here 
achronal will mean spatial or null and causal will mean 
temporal or null. 

A cylinder is said to be trapped, marginal or untrapped 
as V"r is temporal, null or spatial respectively. Assum- 
ing a time orientation, if V*r is future (respectively past) 
causal, then the cylinder is future (respectively past) 
trapped or marginal. On an untrapped or marginal cylin- 
der, an achronal normal direction is outward (respectively 
inward) if r is increasing (respectively decreasing) in that 
direction. This provides a local spatial orientation in an 
untrapped region. A marginal cylinder is outer, degener- 
ate or inner as V 2 r > 0, V 2 r = or V 2 r < respectively. 
A trapping horizon is a hypersurface foliated by marginal 
cylinders. A future (respectively past) outer trapping 
horizon is proposed as the local, dynamical definition of 
a black (respectively white) hole. Trapping horizons de- 
fine what some authors seem to call apparent horizons, 
though the textbook definition of apparent horizon dif- 
fers, being the boundary of the region containing outer 
trapped surfaces in an asymptotically flat spatial hyper- 
surface |@. 
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One may also say that an axial singularity is trapped 
(respectively untrapped) if it has a neighbourhood of 
trapped (respectively untrapped) cylinders, cf. Chiba 
111 . If the singularity is smooth (as a boundary in 
the quotient space) it can alternatively be defined to be 
trapped, marginal or untrapped as V'r is temporal, null 
or spatial respectively. A local formulation of weak cos- 
mic censorship would be that only trapped (or perhaps 
marginal) singularities can form from physically reason- 
able matter and initial data. Strong cosmic censorship 
would require the formation of spatial (or future-null) 
singularities only, where the causal nature of the singu- 
larity is defined by the quotient metric. However, as in 
the spherically symmetric case it is straightforward 
to show that there is an equivalence between these strong 
and weak versions: a smooth axial singularity is spatial, 
null or temporal as it is trapped, marginal or untrapped 
respectively. 

IV. ENERGY 

In Thome's paper on cylindrical energy |?J there is a 
note added in proof to the effect that a certain modi- 
fication of the definition renders it finite in space-time. 
This modified definition, also used by Chiba |ll| , can be 
written as 

E= (l-^ 2 Vr- V tt r)/8 (9) 

with units such that Newton's gravitational constant is 
unity. It transpires that this is indeed the appropriate 
definition of gravitational energy per specific length, for 
reasons developed in the remainder of this article. Princi- 
pally, E is the charge associated with a conserved current 
which combines the energy-momentum of the matter and 
gravitational waves. For the moment, some other prop- 
erties of E will be noted. 

The definition has a similar form to that of the energy 
in spherical symmetry, the Misner-Sharp energy The 
two numerical coefficients are fixed so that E vanishes in 
flat space-time and recovers the mass per specific length 
of a cylinder to leading order in the small-cylinder ap- 
proximation at a regular axis, as shown by Thorne for 
the original definition, the result for E following directly. 
In general, E = 0(p 2 ) at a regular axis, directly from 
the regularity conditions. It also follows that E has the 
correct Newtonian limit, the mass per specific length of 
the cylinder. An example of post-Newtonian behaviour 
is given later for the Levi-Civita metric [O]. 

A cylinder is trapped, marginal or untrapped as E > 
1/8, E = 1/8 or E < 1/8 respectively. This is similar 
to the condition in spherical symmetry but with a 
numerical value arising because E has units of energy 
over length. The value 1 /8 seems to be peculiar to cylin- 
drical symmetry. The corresponding term in E may be 
regarded as a topological term, not present in the corre- 
sponding energy in planar symmetry, which is essentially 
the Hawking energy 



V. TIME 

A canonical flow of time in cylindrical symmetry is 
generated by the vector 

k = [*drf (10) 

where d is the exterior derivative and * the Hodge oper- 
ator of the quotient space, i.e. *d is a curl. The sign of 
the curl can be fixed in an untrapped region so that k is 
future-causal. One may define k equivalently, up to sign, 
by 

/fc-Vr = (11) 
k ■ fc b = -Vr • V s r. (12) 

Then a cylinder is trapped, marginal or untrapped as k 
is spatial, null or temporal respectively. In particular, 
trapping horizons may equivalently be defined as hyper- 
surfaces where k is null. In these and other ways, k is 
analogous to the Killing vector of a stationary space-time 
or the Kodama vector of a spherically symmetric space- 
time @. 

Since the divergence of a vector v normal to the cylin- 
ders is given by V ■ v — r _1 *d*(rv b ), it follows that k is 
covariantly conserved: 

V-/c = 0. (13) 

This conserved current therefore admits a charge 

Q[k] = - I *-k (14) 

which depends only on the boundary <9E of a spatial hy- 
persurface £ with regular axis. Here * is the space-time 
volume 4-form. This charge is found to be the areal vol- 
ume: 

Q[k] = J Vdz. (15) 

The energy-momentum density per specific length of 
the matter, referred to the canonical flow of time, is given 
by the vector 

j[T] = -£- 2 (T-kf (16) 

where T is the energy tensor of the matter in covari- 
ant form. In spherical symmetry, the analogous j[T] is 
also conserved, but this is generally not so in cylindri- 
cal symmetry. The physical reason is that gravitational 
waves carry energy. One might then wonder whether 
there is a definition of energy-momentum for the gravita- 
tional waves, to be added to j[T) such that the combined 
energy-momentum is conserved. It transpires that this is 
so, as shown in the next section. 



3 



VI. ENERGY-MOMENTUM OF 
GRAVITATIONAL WAVES AND THE 
GRAVITATIONAL POTENTIAL 

The energy-momentum density per specific length turns 
out to be 

j=A~ 1 (*dE) t . (17) 

This is the modified version of the vector introduced by 
Thorne 0, also used by Chiba Like k, it is covari- 
antly conserved: 

V • j = 0. (18) 
This conserved current therefore also admits a charge 

Q\j} = - f*-j (19) 

which is found to be the energy: 

Q\j] = J Edz. (20) 

Applying the Einstein equation, j does indeed have a 
term in j [T] , as one can see explicitly in standard coordi- 
nates in the Appendix. Remarkably, the remaining term 
can be written in the same form. That is, two compo- 
nents of the Einstein equation may be written as 

3=j[T]+m) (21) 

where 

j[e] = -r 2 (e-fe)« (22) 

and 

8tt0 = 2Vct)®Vct)~ (V</> • V^ff (23) 

where ® denotes the symmetric tensor product, g is the 
space-time metric and 

4> = - In I. (24) 

Then j[0] is the energy-momentum density per specific 
length of the gravitational waves, which when added to 
j[T] yields the conserved current j. 

One may therefore interpret O as the effective energy 
tensor of the gravitational waves. Such effective energy 
tensors are known in the linearized approximation and 
the short-wave approximation fl4| , but the above result 
is exact: is an invariantly defined tensor. As shown in 
the Appendix, the Einstein, Tolman and Landau-Lifshitz 
energy-momentum complexes as calculated by Rosen and 
Virbhadra Jl5|,[l6| also partially agree with 0, but are not 
invariant tensors. One may recognise as the energy 
tensor of a massless Klein-Gordon field (p. Thus <f> plays 
the role of a potential for the gravitational waves. 



Note that only the components of orthogonal to the 
cylinders contribute to j[0], but it transpires that the 
other components simplify the remaining components of 
the Einstein equation, principally the 2-dimensional wave 
equation 

*d*d 7 = 8tt(T^ + 0£) (25) 

where 7 is defined by the determinant of the quotient 
metric: 

e 27 = ^V^det. (26) 

Although 7 is not invariantly defined, its 2-dimensional 
Laplacian *<i*<i 7 is invariant. Since T£ is the azimuthal 
pressure, the above equation indicates that 0£ acts as 
an azimuthal pressure due to the gravitational waves. 

Another component of the Einstein equation may be 
written as a wave equation for cj> with another invariant 
of the energy tensor T as a source: 

V 2 = 4tt£ (27) 

where 

Q =-tiT-T% + T z z (28) 

and the trace refers to the quotient space. This may also 
be written with T + replacing T, since 

trO = (29) 
= 9*. (30) 

The same is true for the one component of the Einstein 
equation yet to be given. Thus is an effective energy 
tensor in a general sense: writing the Einstein equation 
with T + on one side, instead of the usual T, simplifies 
the other side. 

The wave equation for <f> is recognisable as a relativis- 
tic version of the Poisson equation of Newtonian gravity, 
with <fi playing the role of the Newtonian gravitational 
potential. This was pointed out by Melvin p||, though 
with the opposite sign for the potential. The sign is de- 
termined by the fact that g reduces to the density in the 
Newtonian limit. In summary, the gravitational potential 
(j) both generalizes the Newtonian potential and acts as 
a potential for the gravitational waves. 

VII. UNIFIED FIRST LAW 

The first law is the energy-balance equation expressing 
the gradient of E according to the Einstein equation. 
This can be written in a similar form to that of spherical 
symmetry, 

dE = Aip + r 2 wdV (31) 
in terms of the work density (an energy density) 
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w = -trT/2 
and the energy flux per specific length 
tP = i/j[T} + ip[Q] 



(32) 



(33) 



which has been divided into contributions from the mat- 
ter and the gravitational waves: 



f/>[Q] 



wVr) 



(34) 
(35) 



These energy fluxes are essentially duals of the energy- 
momentum densities: 



f[T] 



*iP[T] - 
*ip[0]. 



(36) 
(37) 



The terms in the first law involving ip and w may be in- 
terpreted as energy-supply and work terms respectively, 
analogous to the heat supply and work in the classical 
first law of thermodynamics p7| . This particular divi- 
sion into two terms is motivated by properties of ip, de- 
scribed below, which make it analogous to heat flux in 
thermodynamics. The unified first law was so called be- 
cause it was shown in spherical symmetry that projecting 
it along the flow of a thermodynamic fluid yields a first 
law of relativistic thermodynamics, while projecting it 
along a trapping horizon yields a first law of black-hole 
dynamics ||. In cylindrical symmetry, the above first 
law is also unified in the sense that it includes the energy 
flux of the gravitational waves. It may also be applied to 
cosmic strings, as shown later. 

The energy density and fluxes have similar properties 
to those of the spherically symmetric case ||. Firstly, 
ip$[@] is past (respectively future) causal in future (re- 
spectively past) trapped regions, and outward achronal in 
untrapped regions. This last property of ip$ [&] is analo- 
gous to that of heat flux, which is a spatial vector in ther- 
modynamics |rij| . The same causal properties hold for 
ip$ [T] assuming the null energy condition. Then £ • ip > 
where £ is an outward achronal vector. Also w > as- 
suming the dominant energy condition. Then the first 
law implies 



C-V£>0. 



(38) 



Thus E is non-decreasing in any outward achronal direc- 
tion in an untrapped region. This monotonicity property 
immediately gives a positivity property: on an untrapped 
achronal hypersurface with regular axis, 



E > 0. 



(39) 



Similarly, on an untrapped achronal hypersurface with 
a cosmic string of energy per specific length e, defined 
later, 



This type of argument can be used to give several inequal- 
ities for black holes ifTsjf . If the space-time is asymptot- 
ically flat at spatial or null infinity, one can define the 
asymptotic mass per specific length to be the limit of 
E, if it exists. Then the same inequalities hold for the 
asymptotic mass. The monotonicity property reduces at 
null infinity to a mass-loss property, with ip reducing to 
a Bondi-like flux per specific length . Generally, ip is 
the outward flux minus the inward flux. 



VIII. SURFACE GRAVITY 



A dynamic surface graviti 
spherically symmetric case J 



k may be defined as in the 
I by 



k - (V Afc b ) =&Vr 



(41) 



where A denotes the antisymmetric tensor product. That 
is, the 1-forms on each side of the equation are propor- 
tional and k is defined as the proportionality constant. 
Since k — ±V"r on a trapping horizon, this reduces on 
such a horizon to 



k ■ (V A k b ) = ±ink b 



(42) 



which is analogous to the usual definition of stationary 
surface gravity, with k replacing the stationary Killing 
vector. Then 



2£k = *dk — *d*dr. 



(43) 



Since V 2 



*d*dr on a marginal surface, this shows 



that a trapping horizon is outer, degenerate or inner as 
K>0,K = 0or«<0 respectively. This confirms a 
desired property of surface gravity, that it should vanish 
for degenerate black holes. As in spherical symmetry ||, 
k is defined everywhere in the space-time, not just on 
horizons. 

The remaining component of the Einstein equation can 
be written as a 2-dimensional wave equation for r: 



Then 



vd*dr — —\Qitrw. 



Sirpw. 



(44) 



(45) 



Therefore black or white holes, as defined by outer trap- 
ping horizons, are not consistent with the dominant en- 
ergy condition. This is implied by the general topology 
law for black holes Q| and reflects an unphysical aspect of 
cylindrical black holes. However, note that it is possible 
to maintain the null energy condition while violating the 
dominant energy condition. The simplest example is a 
negative cosmological constant A < 0: 8irT = —Kg. 



E > e. 



(40) 
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IX. LAWS OF BLACK-HOLE DYNAMICS 



X. STATIC EXAMPLES 



Projecting the first law along a trapping horizon yields 

= £E> = — +r 1 wV (46) 

where /' = £ • V/ in terms of a vector £ tangent to the 
horizon; the quickest way to derive this uses the previ- 
ous displayed equation. Corresponding equations in the 
spherically symmetric case [|| and the general case |2p| l 
have been called the first law of black-hole dynamics, as 
they take the same form as the familar first law of black- 
hole statics, which actually is more like a Gibbs equation 
if one takes the thermodynamic parallel seriously §IJ§ . In 
any case, the surface gravity and area occur as expected. 
By analogy to the original definition of entropy by Clau- 
sius [^lj, this identifies A as a specific entropy and k 
as a temperature, up to factors. The fact that the en- 
tropy arises from the energy-supply term again indicates 
that ij) is analogous to heat flux, which is proportional to 
entropy flux in conventional thermodynamics, effectively 
defining temperature |T7j ]. 

The second law of black-hole dynamics jl| reads 

A! > (47) 

assuming the null energy condition, where the outward 
orientation of the future outer trapping horizon is such 
that the outward null expansion vanishes. A trivial ze- 
roth law expresses the constancy of k over a cylinder. Of 
course, k generally varies along a trapping horizon, unlike 
the static case. The signature law of black-hole dynam- 
ics 0| states that an outer trapping horizon is achronal, 
again assuming the null energy condition. Without the 
null energy condition, outer trapping horizons may also 
be used to define traversible wormholes 0]. Thus one 
may also investigate the interaction of wormholes and 
gravitational waves. 

As indicated above, black holes in cylindrical sym- 
metry have some peculiar features. Firstly, their exis- 
tence requires violation of the dominant energy condi- 
tion, whereas their usual properties require the null en- 
ergy condition. In axial symmetry, the dominant energy 
condition therefore implies a limit on the formation of 
arbitrarily spindly black holes, as also argued by Thorne 
jl4],^7| ■ Secondly, a black hole has constant energy per 
specific length, E = 1/8. In particular, there is an en- 
ergy gap between black holes and flat space-time. Thus 
it is impossible to smoothly develop a trapped region 
around an initially regular axis, where E = 0. Unless 
E is badly behaved, e.g. jumping or becoming infinite, 
gravitational collapse to a black hole first requires the 
initially regular axis to develop into an axial singularity 
with < E < 1/8. This is essentially a cosmic string, as 
the following examples indicate. 



Static metrics can be locally specified by (£, E, F) as 
functions of r, where F = — £t • ^ in terms of the static 
Killing vector £ t . Then the line-element is 

ds 2 = £ 2 dz 2 + r 2 r 2 dip 2 + r 2 (l - 8E)- 1 dr 2 - Fdt 2 . 

(48) 

A simple example is the metric of a static cosmic string, 
which may be written as 

ds 2 = dz 2 + r 2 d(p 2 + (1 - 8s)- 1 {dr 2 - dt 2 ) (49) 

where the constant e is interpreted as the energy per 
specific length of the string, since E = e. The space- 
time is locally flat, but has a conical axial singularity 
unless e vanishes. Normally the case < e < 1/8 is 
considered, so that there is a positive angular deficit 9 = 
27r(l-(l-8e) 1 / 2 ). For small e, 9 = 8ne + 0(e 2 ), agreeing 
to lowest order with the usual expression for the mass or 
energy of a cosmic string. One may also consider e < 
and e > 1/8, the latter being a trapped, spatial string. 

Cosmic strings are weak axial singularities in that E 
is finite. An example of a stron g ax ial singularity is pro- 
vided by the Levi-Civita metric |12| . This is the relativis- 
tic analogue of the Newtonian gravitational field outside 
an infinite cylinder, given by the Newtonian gravitational 
potential 2mlnr, where m is the mass per specific length. 
In Einstein gravity, the unique static vacuum solution has 
two constants, due to the additional possibility of a su- 
perposed cosmic string. The solution may be written, cf. 
Thorne 0, as 

ds 2 = r~ im dz 2 + r im+2 d V 2 

+ (1 - 8m - 8£ )-i r 4m ( 1+2m ) (dr 2 - dt 2 ) (50) 

where the constants have been arranged so that m is the 
mass per specific length and e is the energy per specific 
length. That is, if one expands 

E = 1/8 - (1/8 - m - e)r- 8m2 (51) 

by assuming such units in the Newtonian limit, i.e. 
(r, to, e, E) i — ► (r, c~ 2 m, c~ 4 e, c^^E) where c is the speed 
of light, one finds 

E = mc 2 +e + m 2 lnr + 0{cr 2 ). (52) 

The third term has the same form as the expression 
TO0/2 for Newtonian gravitational self-energy per spe- 
cific length, where <fi = 2m In r as in the Newtonian case. 
Thus E correctly encodes the mass of the line source, 
its gravitational self-energy and the energy of the super- 
posed string, which may also be regarded as analogous to 
internal energy. This allows a curious quasi-Newtonian 
distinction between energy and mass, only the latter pro- 
ducing a non-flat gravitational field. 
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Finally, a simple example illustrating strings and black 
holes is the static metric given by 

ds 2 = dz 2 + r 2 dif 2 + (1 - 8e + aV) -1 ^ 2 

-(l-8e + a 2 r 2 )dt 2 . (53) 

This is the trivial lift of the spinless 3-dimensional BTZ 
black hole [^3| , which is a solution for a negative cosmo- 
logical constant A = —a 2 , reducing to the 3-dimensional 
anti de Sitter solution if e — 0. Additional pressure terms 
are required in the 4-dimensional case p4| ]. If a = 0, a 
static cosmic string is recovered, as above with a different 
scaling of t. Generally 



E = e-a 2 r 2 , 



(54) 



so that there is a string of energy e per specific length, 
with the other term in E being the contribution of 
the cosmological constant, analogously to the spherically 
symmetric case 25 . If e < 1/8, the string is naked. If 
e > 1/8 there is a trapping horizon at r — a _1 (8e — l) 1 / 2 . 
This is an outer trapping horizon with surface gravity 
k = a 2 r. Thus there is a black hole with a similar global 
structure to the Schwarzschild black hole. In particular, 
the string inside the black hole is spatial and trapped. 

It was originally argued that the Euclidean version of 
the BTZ black hole yields a temperature k/2it, as ex- 
pected, and an entropy Airr [^3| . However, this involved a 
mass defined as 8e — 1. Repeating the argument with e as 
the mass yields an entropy 7rr/2, or in the 4-dimensional 
case a specific entropy A/ A. This expected value for en- 
tropy has since been confirmed by statistical-mechanical 
methods [Eft. 



XI. COSMIC STRINGS 

Since cosmic strings may be a feature in gravitational 
collapse in cylindrical symmetry, it may be useful to con- 
sider dynamic strings. Here a cosmic string will be de- 
fined as a smooth axial singularity for which E has a 
finite limit, unique at each point (in the boundary of the 
quotient space). This limit 



£ = lim E 

r— >o 



(55) 



is the energy per specific length of the string. Then a 
string is trapped and spatial if £ > 1/8, marginal and 
null if £ = 1/8, and untrapped and temporal if e < 1/8, 
cf. Chiba flrifl. In particular, unlike the spherically sym- 
metric case there is now a range of positive energy, 
< £ < 1/8, for which axial singularities are locally 
naked. Thus if gravitational collapse to a black hole is to 
occur, then provided £ is smooth, the initially regular axis 
must first develop into a cosmic string with small energy 
£. Then only if e increases past 1/8 will a trapped region 
form, hiding the string. Until or unless that happens, the 
string is a locally naked singularity. 



This reflects the ease with which naked sing ularities 
seem to form in cylindrical symmetry [Tl|j27| , p8|| or axial 
symmetry [|l^,|9[|o), at least for dust. This energy gap 
between flat space-time and black holes can be seen as the 
physical basis for the apparent failure of cosmic censor- 
ship in cylindrical symmetry. However, this could also be 
ascribed to the pathological nature of dust, since similar 
studies for a fluid with pressure indicate that a bounce 
occurs instead p7| , pl| . Another possible resolution sug- 
gested by Nakamura et al. |}2| is that the energy of the 
incipient singularity is dispersed by gravitational waves, 
presumably preventing its formation or minimizing its 
effect. Studies of gravitational- wave emission prior to 
singularity formation have yielded quite different results: 
negligible emission |ll| , ^9| -pll or a strong burst [^8) . 

These issues can be studied in the context of cosmic 
strings. The first law projected along a string reads 



£ = lim Ak 



(56) 



where f = k- V/. This is the first law for cosmic strings, 
relating the increase or decrease of £ with the absorp- 
tion or emission, respectively, of gravitational waves or 
matter. Clearly e is constant if the flux k ■ tp is finite 
at the string, as occurs for Einstein-Rosen waves and a 
massless Klein- Gordon field. Such radiation seems not to 
feel the presence of the string at all, the incoming waves 
passing through it and re-emerging, with no gravitational 
collapse. On the other hand, if gravitational collapse to 
an axial singularity does occur, this is a cosmic string 
with time-dependent energy if the flux k ■ ifi is 0(r _1 ), 
i.e. for finite energy supply Ak ■ tp. A stronger axial sin- 
gularity could conceivably occur, with effectively infinite 
energy. Either way, predictability is lost unless there is a 
prescription for the emitted radiation. 

One may take the view that naked singularities are not 
problematic if predictability can be maintained. For in- 
stance, Clarke |33j has defined generalized hyperbolicity 
in terms of existence and uniqueness of solutions to the 
wave equation, and Vickers & Wilson [Q have shown 
that a static cosmic string is hyperbolic in such a sense. 
To handle a dynamic cosmic string, however, it seems 
necessary to impose a boundary condition at the string. 

Without an internal model for the string, the simplest 
prescription is the reflecting boundary condition, 



lim Ak ■ ib\ 







(57) 



applied separately for O and any types of matter con- 
tributing to T . (If only the combined flux ip were so 
constrained, energy transfer between gravitational waves 
and matter would be allowed at the string, again remov- 
ing predictability). This boundary condition means that 
any radiation striking the string is immediately reflected, 
leaving the energy of the string constant: 



0. 



(58) 
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This is automatically obeyed for Einstein-Rosen waves, 
but would be inconsistent with gravitational collapse 
as outlined above; some absorption of radiation by the 
string is necessary for its energy to increase from the zero 
of initial regularity. An alternative boundary condition 
is the absorbing condition 



lim Ak ib- 







(59) 



where <9_ is the future-inward null direction, as in the 
Appendix. In this case 



e > 



(60) 



assuming the null energy condition, with the string ab- 
sorbing all matter and radiation which strike it and emit- 
ting nothing. One might think of such irreversibility as 
a second law for cosmic strings. 

Another possibility is a linear combination of the above 
conditions: 



lim A(ek+ip+ + k~ip-)\. ..1=0 

r— h0 



(61) 



where d+ is the future-outward null direction and the 
reflection coefficient e is specified between and 1, for 
no reflection and total reflection respectively. There are 
yet more possibilities, for instance in a dust collapse one 
might impose an absorbing condition for the dust, but 
a reflecting condition for the gravitational waves. The 
question of whether such naked singularities are strong 
sources of gravitational waves is presumably dependent 
on the choice of boundary condition. 



XII. CONCLUSION 

In cylindrically symmetric Einstein gravity, the phys- 
ical quantities may be summarized as follows, omitting 
the recurring qualification per specific length. The kine- 
matical or gravitational quantities are the specific area 
A (or r = A/2ir), the canonical time vector k, the grav- 
itational potential <p (or the specific length t = e~^, or 
the circumferential radius p — r/£), the energy E and the 
surface gravity n. The effective energy tensor O is defined 
in terms of <f>, leading also to the energy-momentum j [O] 
(or the energy flux ^[©D of the gravitational waves. The 
matter quantities are the pressures and T£, the en- 
ergy density w (or g) and the energy- momentum j [T] (or 
the energy flux V'PDi au °f which are invariants of the 
energy tensor T. 

The five independent components of the Einstein equa- 
tion may be written as the first law plus three wave equa- 
tions, each with sources given by invariants of T + 0. 
They are collected here as 



dE = AtP + r 2 wdV 
*d*dr = — Wirrw 
*d*dj 
V 2 (j) = 4ttq 



(62) 
(63) 
(64) 

(65) 



Instead of the first law, one may use the equation ex- 
pressing j, the energy-momentum associated with E, as 
j [T + O] . Energy conservation is expressed by the covari- 
ant conservation law V • j = 0. In summary, the Einstein 
equation takes an invariant form, with the combined en- 
ergy tensor T + as a source. In this sense, is an 
effective energy tensor for the gravitational waves. 

At a time when detection of gravitational waves is ex- 
pected to be imminent, it may seem odd to recall that 
there was a historical debate about whether gravitational 
waves carry energy, or are well defined at all in Einstein 
theory. Einstein-Rosen waves played a significant role in 
this debate 15 lr^ , |35| , |36[ | . There were several attempts to 
construct pseudotensors for the energy, which are usu- 
ally rejected nowadays due to their non-invariant nature. 
The modern viewpoint seems to be that expressed by 
Misner, Thorne & Wheeler Jl4|], that gravitational waves 
and their energy may be defined only in an approximate 
or averaged sense. Thus it may come as a surprise that a 
genuine tensor exists which encodes the energy of the 
gravitational waves in cylindrical symmetry. One may 
partly ascribe this to the symmetry, since averaging is 
trivial over the cylinders. However, there has been no 
need to average over the remaining dimensions, or to as- 
sume small amplitude or wavelength, as can be seen ex- 
plicitly for Einstein-Rosen waves in the Appendix. The 
discovery of should at least resolve the debate about 
the energy of Einstein-Rosen waves. 

Basic laws of black-hole dynamics have been given, 
including a first law with the expected term involving 
area and surface gravity. There is already evidence that 
these determine entropy and temperature for static black 
holes p6| ) and it may be conjectured that this general- 
izes to dynamic black holes. Recently there has been 
considerable interest in anothergeneralization of static 
black holes, isolated horizons pT\ , which are a certain 
type of null trapping horizon. In particular, a quantum- 
geometrical derivation of black-hole entropy can be given 
for such local horizons (3^] . 

Finally, there is a formula for the energy which unifies 
cylindrical and spherical symmetry: 



M 



C_ 

An 



1 - 



VA • V>U 



Kir- ) 



(66) 



In spherical symmetry M is the Misner-Sharp energy 
if A is the area and C the circumference of the spheres. 
In cylindrical symmetry M is the specific energy IE if 
r = nl/2. In each case, a sphere or cylinder is trapped 
if and only if 



r < AwM. 



(67) 



This is exactly the inequality in the hoop conjecture of 
Thorne 27|. It is to be hoped that the theory developed 
here can be generalized to axial symmetry, so as to have 
a fully dynamical description of rotating black holes and 
co-rotating black-hole collisions. 
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XIII. APPENDIX 

In cylindrical symmetry as defined here, the line- 
element locally takes the form 



ds 2 



tdz 2 + p 2 dip 2 + hijdx'dx 1 



(68) 



where x l are coordinates orthogonal to the cylinders and 
the metric functions (p,£,h) are independent of ((p,z). 
Using double-null coordinates x ± and definitions in the 
text, this takes the form 



ds 2 



~^d 



_ , 2d, 2 

z + e v r 



d^ 2 - 2e 2 ^dx+dx~ 



(69) 



The remaining coordinate freedom consists of the null 
diffcomorphisms x ± i— ► ^(ar^) and constant linear trans- 
formations z i ^ Bz + C, with ip identified at and 2n. 
Then for a function / and a 1-form a, 



(df)± = d±f 
(*a)± = T«± 



*da 



-27-20 



(d-a + — 9 + a_) 
V A a = da/2 

V 2 / = *d*df + r^V'r ■ V/ 



(70) 
(71) 
(72) 
(73) 
(74) 



= -e- 2 ^ ( 2d + d-f + O+rB-f + B-ra+r ^ 



where the orientation of * has been locally fixed so that 
d + is future-outward and 9_ is future-inward on un- 
trapped or marginal cylinders. Then various quantities 
defined in the text may be written as 



E= {l + 2e- 2 ^d+rd-r)/8 



J 



-1-20-27 



(76) 
(77) 

d T E (78) 
j±lT] = ±e- 2 ^(T TT d ± r - T+_9 T r) (79) 
j±[e] ^±e- 24, -^(d T cb) 2 d ± r/ATT (80) 

e ±± = {d ± <t>) 2 /AK (si) 

(82) 
(83) 
(84) 
(85) 
(86) 
(87) 



e+_ = o 

* = e^+^rdz A dtp A dx + A dx~ 
w = e- 2 ^ 27 T+_ 
i>±[T\ = -e- 2 ~*T ± ±d T r 
^±[6] = -e- 2l {d±4>) 2 d T r/An 
k = -e-^-^d+d.r. 

The non-zero Christoffel symbols are 



r* = 

zz 



e" 27 r(a T r + rd T <j)) 



(89) 
(90) 



r| ± = 2(9±0 + a± 7 ) 
r± =2{d^4> + d Tl ) 



r 



Z± 



(91) 
(92) 
(93) 
(94) 



and the non-zero components of the Ricci tensor are 

R zz = -e-^- 2 ~<x 

(2d + d-<j> + r- 1 (d + rd-(p + d-rd + (j))) (95) 

{2d+d-r + 2rd+d-<i) + d + rd-<i) + d-rd + <t)) (96) 

R±± = r- 1 (2d±rd±j - d±d±r) - 2{d±4>) 2 (97) 

R+- = -2d+d-<f) - 2<9 + <9_7 - 2d + <pd-(j) 

-r-^d+d-r + d+rd-^ + d-rd+cf)). (98) 

The Einstein equation may then be written as 

d±d±r - 2d ± rd ±1 + 2r(<9±</>) 2 = -8irrT±± (99) 
d+d-r = 8nrT+- (100) 
9+5-7 + d+4>d-4> = -Awe^+^Tg (101) 

2d+d-4> + r- 1 (d+rd„(f> + d-rd+cj)) 

= 47r(e 20+27 (r^-T;)-2T + _). (102) 

It follows that 



d ± E = 2nre~ 2 ~> ( T + _d±r - ( T±± + ) 5 T r 



4tt 



(103) 



The null energy condition requires T±± > and the dom- 
inant energy condition requires T_| > 0. The equations 

in the main text follow straightforwardly from the above 
expressions. 
Writing 



V2x ± = x°±x 1 



(104) 



puts the metric in the form of Melvin |J and Thorne [0 . 
The Einstein equation in these coordinates is found by 
linear combinations to be 



r-r" = -87re 27+2 ^r(T + I?) 



" j.'2 j.2 



57re 2 7 +20 T ^ 



(105) 
(106) 

4> - (j)" + r- l {r4> - r'cj)') = Aire^+^x 

(T ° + Tl + T£ - 1?) (107) 

(r' 2 - f 2 ) 7 ' = ~8^e 27+2 ^(r'T ° + fTj) 

+rr'(<j) 2 + 4>' 2 ) - 2rr<j)(j)' + r'r" - ff (108) 
(r' 2 - f 2 ) 7 = 87re 2 T+ 2 *r(r'r 1 - rT^) 

-rr(4> 2 + <f)' 2 ) + 2rr>0' - rr + r'r' (109) 

where / = df/dx° and /' = df/dx 1 in this appendix 
only. This is essentially the form given by Melvin jf| and 
Thorne except for misprints in their versions of the 
last equation. 
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If w = 0, e.g. in vacuo, one may solve the wave equa- 
tion for r by taking x 1 — r locally in an untrapped region, 
in which case the metric takes the Einstein-Rosen form 
j|. Note however that what most authors call i> has 
been replaced with —<f>, since <j) generalizes the Newto- 
nian gravitational potential. The energy tensor of the 
gravitational waves takes the form 

87t6 = (0 2 + (j> 12 ) (dt (g>dt + dr(g> dr) + 2^'dt ® dr 

+ {<j) 2 ~ (j)' 2 )e- 2l {e~ H dz ®dz + r 2 dip <g> dip) (110) 

where t = x . The energy density Ooo and energy flux 
Ooi agree up to factors with the Einstein, Tolman and 
Landau-Lifshitz energy-momentum complexes as calcu- 
lated by Rosen and Virbhadra |Hj|n|. The pressure 
terms seem to be new. 

Einstein-Rosen waves are given by the Fourier-Bessel 
modes 



A{u))e Mt J {u;r)dLu 



(111) 



where Jo is the zeroth-order Bessel function and the func- 
tion A is complex. This solves the wave equation for 
<j> in vacuo, the complete solution to the vacuum Ein- 
stein equation following by quadrature for 7. For a single 
mode, substitution shows that depends on the squares 
of the amplitude \ A\ and frequency uj. In particular, the 
energy density Ooo at the axis oscillates between and 
|^4| 2 w 2 /87r, as for an electromagnetic wave. Unlike the 
linearized or shortwave approximations [ p"4| |, \A\ and u> 
can be arbitrarily large. These fully non-linear waves 
still admit an energy tensor O. 
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